"Chain scenario" for Josephson tunneling with vr-shift in YBa2Cu307 
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We point out that all current Josephson-junction experiments 
probing directly the symmetry of the superconducting state 
in YBa2Cu307, can be interpreted in terms of the bilayer 
antiferromagnetic spin fluctuation model, which renders the 
superconducting state with the order parameters of extended 
s symmetry, but with the opposite signs in the bonding and 
antibonding Cu-O plane bands. The essential part of our in- 
terpretation includes the Cu-O chain band which would have 
the order parameter of the same sign as antibonding plane 
band. We show that in this case net Josephson currents along 
and perpendicular to the chains have the phase shift equal to 



In the last year, starting with the pioneering work of 
Wollman ct al substantial progress has been made 
|^,^|^J^,^,^ in probing the symmetry of superconducting 
state in YBa2Cu307 by mean of Josephson tunneling. In 
all these experiments except of |Q] the relative phase of 
the tunneling currents in YBCO contacts parallel to a 
and to b crystallographic axes was measured. In most 
cases it was found that the phases are opposite, as ex- 
pected for instance for d^2_y2. To the contrary, in Ref. 
Q tunneling current parallel to c was measured, which 
for pure dj;2_y2 is expected to vanish and non-zero, 
although small, value was found. 

Interpretation of the experiments is addi- 

tionally obscured by the fact that the only object studied 
so far has been orthorhombic YBa2Cu307, where a. d + s 
state is formally allowed and one can speak only about 
the weight of d— or s— components. Some authors 
suggested that a strong anisotropy of the Fermi surface 
can explain the edge-contact experiments even without 
large d— component. However, the underlying assump- 
tion is that the plane electrons themselves are subject to 
strong orthorhombic effects, while both in calculations 
pO| and in the experiment [ pl] | the main manifestation 
of orthorhombicity is the presence of the chains, while 
the planes themselves remain fairly tetragonal. This fact 
cannot be neglected when judging about pairing symme- 
try (see, e.g., Ref. 

In a previous work |p^ , we noticed that if the order 
parameters (OP) in chain and plane bands had oppo- 
site sign and if the tunneling current along the chains 
were dominated by the chain band, this could explain 
the Josephson experiments in YBCO (this suggestion has 
been recently elaborated on by others [0). While in 



Ref. p^ ] a number of reasons have been proposed for the 
sign reversal of the OP, none of Refs. p3| , p^ suggested 
any microscopical reason for the tunneling current being 
dominated by chains. 

In this Letter we propose another, quantitative "chain 
scenario" for the above-mentioned Josephson tunneling 
experiments. This scenario is based on a recently pro- 
posed bilayer antiferromagnetic spin fluctuation model 
for superconductivity in YBCO [l5|] where the symmetry 
of the pairing state for the plane electrons is such that 
the bonding and antibonding plane bands have OP of the 
opposite signs while angular symmetry is extended s. We 
will argue that within this model, if the detailed Fermi 
surface properties of YBa2Cu307, as calculated in LDA 
and measured experimentally, are taken into account, it 
follows that the net tunneling currents along a and along 
b must have opposite signs. 

The Fermi surface of YBa2Cu307 is believed to con- 
sist of four sheets: Two plane bands, which are bonding 
(b) and antibonding (a) combinations of the individual 
planes' states (in the calculations they are well splitted 
by the energy which ranges from 0.05 to 0.5 eV, depend- 
ing on the wave vector), the chain (c) band, and a small 
pocket formed mainly by apical oxygen states (which is 
not discussed here). The Fermi surface of YBa2Cu307, 
as calculated by Andersen et al, jl^ is shown in Fig. 
1. The c band has been detected by positron annihila- 
tion technique and the corresponding sheet of the Fermi 
surface is in perfect agreement with LDA calculations. 
According to calculations [^, this band is very light, 
so that its contribution to the total density of states is 
small (~ 15%), while its contribution in the plasma fre- 
quency Wp cx N{0)vp is considerable (~ 50%). These 



1 



finding are confirmed by the experiment: Maximal Fermi 
velocity was calculated to be ~ 6 x 10'' cm/s and 
corresponds to the point where the chain Fermi surface 
crosses the F— Y line. This value agrees well with the 
Raman experiments JlSt . Calculated plasma frequency 
anisotropy uj^y/uj^^ ~ 1.75, as discussed in Ref. is in 
agreement with the optical and transport measurements. 

Band a is, according to calculations, rather heavy. In 
particular, it has extended van Hove singularities which 
have also been discovered experimentally |l^j2^ . Band h 
is light again. Both a and b bands are nearly tetragonal. 
Their relative contribution to the normal-state transport 
is defined by the partial plasma frequencies (Table I). 
Importantly, bands a and c at = can cross by sym- 
metry, for instance they are degenerate with e — Ep at 
q=(« 0.87r/a,R:! 0.27r/6,0). For all 7^ these bands 
hybridize. This is the reason for YBCO being the most 
three-dimensional of all high- Tc cuprates. An extremal 
orbit in = plane, which appears because of the a — c 
hybridization, has been seen in de Haas-van Alphen ex- 
periments |2l[| . These facts provide indirect support for 
the calculations as regards a — b splitting and a — c hy- 
bridization. 

Now we make link to the above mentioned model for 
the superconducting state, suggested in Ref. The 
key feature of the model is that the bands of the differ- 
ent parity, a and b, have OP of the opposite signs. The 
sign of the OP in band c was not discussed in Ref. [ p^ . 
Apparently, because this band hybridizes with band a, 
but not band b, one can assume that c and a have OP of 
the same sign, while b has OP of the opposite sign. How 
can this fact manifest itself in Josephson tunneling? 

To answer this question, let us consider the tunneling 
currents between two superconductors, L and R, each 
having several conducting bands, labeled by subscripts 
For simplicity let us assume that the OP A.^ in in- 
dividulal bands are isotropic. In order to evaluate the 
Josephson current between each of these bands we make 
use of a standard formalism of the integrated over en- 
ergy quasiclassical Green functions which should be 
completed by the boundary conditions on both sides 
of the barrier. Assuming that the tunnel junction trans- 
parency is small -D(p) <C 1 one can proceed perturba- 
tively and expand the boundary conditions ]2^ ] in powers 
of D. Keeping only the linear terms one gets 

gu{p) - 9Li{-p) = ^Aj(/j?j(p)/Li(p)-/Li(p)/flj(p)), 

where the functions g and / are respectively the nor- 
mal and anomalous quasiclassical Green functions on the 
left- or on the righthand side of the barrier. As we 
are interested only in the linear in D(p) contribution 
to the Josephson current it is sufficient to substitute for 
the anomalous Green functions the unperturbed values. 



fLi{B.j) = ^Li{Rf)/^/\^Li{Bj)\'^ +^ni- Then, using 

J = -i2TT)-^ieTY, f {d^S/vF)^Fip)9{p,0Jm), 

we arrive at the general expression for the Josephson cur- 
rent between the bands i and j {Jtot = "^ij)'- 
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Here Rij is the normal state resistances of a tunnel 
junction for the bands (i, j), and is the maximal of the 
two resistances, Rl and Rj^ : 
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v„ is the projection of the Fermi velocity vp on the di- 
rection normal to the junction plane, dS is an element 
of the Fermi surface for the corresponding band. Eq. 
|l| is a straightforward generalization of the well known 
Ambegaokar-Baratoff result to the case of several 
conducting bands. 

Further simplification of Eq. |l]. takes place at low 
temperatures T <^Tc 
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A, log4|A,/A,|/ei?,,, 
7r(A,Aj)/e%(|A,| + |A,|), 
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(A, < A,) 
(A, « A,) 



where Kit) is the complete elliptic integral. The first 
case corresponds to a contact YBa2Cu307— conventional 
superconductor, and the second to a grain boundary con- 
tact. 

Estimate of the effective barrier transparencies Dij in 
Eq.^ can be obtained for certain models of the potential 
barrier U {x) between L and R. Let us consider the case of 
a specular reflecting barrier U{x) = Uo5{x — xq)- Trans- 
mission probability for a quasiparticle from the band i in 
L to tunnel into the band j in R is found by matching 
solutions of Schrodingcr equations on both sides of the 
barrier using boundary conditions at the barrier pq] 
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The second condition is conservation of the probability 
current J{x) — ~ilm{'$*d'i>/dx) /2m{x). It is important 
to note that, as was shown in Refs. rriij are effec- 
tive band masses of quasiparticlcs in L and R , and are 
neither the bare electron mass mo , nor the masses renor- 
malizcd by many-body correlation effects, i.e., essentially 
the LDA band masses. As a result, effective barrier trans- 
parency coefficient Dij in Eq.^ is determined by the band 
velocities: 
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In the low transparency limit, Uq 3> we have Dij — 
,nVRj^m where Dq is a constant. Thus, for a con- 
ventional superconductor on the righthand side, we have 
Ar < Aa,b,c, and: 



Ja • Jb • Jc 
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From Eq.^ we observe that band c does not contribute 
into the tunneling current in the x direction (i.e., per- 
pendicular to the chains), which is quite natural. Corre- 
spondingly, the total current along x is 



Jx^Ja + Jb = \Ja\ - \Jb\ 

Jy^Ja + Jb + Jc= \ Ja\ - \Jb\ + I Jc 



(7) 



. Substituting the values for v from Table 1 in Eq. |g we 
get 



Ja : Jfc : Jc ~ 1 : 2 : 2.2 



(8) 



Now we observe that \Ja\ < |Jb|, while |Jc + Ja| > |Jb|, 
unless \Aa\, \Ab\, and |Ac| differ drastically (so that the 
log terms in Eq.^ become important) . To check this pos- 
sibility, let us come back to the bilayer antiferromagnetic 
spin fluctuation model of Ref. The essence of the 

model is that the coupling interaction is interband a — b 
interaction. A non-essential feature of the model was that 
the densities of states in both bands were assumed equal. 
To estimate the effect of Na being about twice larger that 
Nb, let us consider the weak coupling limit for a BCS su- 
perconductor near Tc with interband interaction only. In 
this limit, 

Aa = const ■ VabNbAb Ab = const ■ VabNaAa, (9) 

where V is the pairing interaction, and we obtain for the 
ratio of the gaps |Aa/A;,| = y^Nb/Na, i.e., counterintu- 
itively, the band with the smaller density of states, in our 
case, bonding band, develops a larger gap. Including Ac 
in Eqs.p assuming Nc <C Na, Vbc ~ Kc ~ 0, we obtain 

\AJAb\ = Vac/Vab- 

Thus, we can safely exclude the possibihty of |Af,| be- 
ing too small, but there are no arguments within the 
chosen model that Ac cannot be arbitrary small. How- 
ever, there are various experimental indications that this 
is not the case, for instance optical experiments of Bauer 
et al p6t , who compared normal/superconducting opti- 
cal conductivity ratios for pure YBCO and YBCO doped 
with Fe (which substitute Cu in the chains). They found 
that the main effect of doping was that a gaplike struc- 
tures at a; ~ 150 cm^^ shifts down to ~ 50 cm~^, while 



the maximal gap at w ~ 300 cm~^ does not change. The 
lower energy can be naturally interpreted as the chain 
gap and the higher one as the plane gap. Thus, one can 
be confident that relations between J^, Jb, and Jc indeed 
hold. 

Now, since Jc and J a are in-phase, and Jb is out-of- 
phase with both of them, it becomes obvious that the 
net current along x and along y should have a phase 
shift TT. Such a state is indistinguishable from d^'^-y'^ 
state for those experiments which probe phase difference 
for two edge contacts; however, for the tunnel current 
perpendicular to the planes the model correctly gives a 
non-zero value. 

The discussion above is relevant for the experiments 
like Refs. which deal with YBCO-conventional 

superconductor contacts. Let us now discuss the case of 
YBCO- YBCO contacts ||] (grain-boundary junctions). 
The condition for a 7r-contact is now | JacI ~ | Jhc| + 
I Jaa| + I JbbI - 2| Jah| < 0, aualogous to Eq||. Accord- 
ing to Eqs. |,|, the corresponding currents are Jij oc 
ViVjAiAj/{\Ai\ + \Aj\). To estimate currents, let us use 
the relations between the order parameters: Aa/ Ab = 
y/Nb/Na and introduce |Ac/Ab| 
the condition is 

aVaVcv^ avbVc , v? , 
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a. Then 



Nb 2VaVbVNb 
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< 0. 



Substituting data from Table I, Va '■ Vb '■ Vc ^ ^ '■ 2 : 2, 
Na : Nb : Nc ^ 2 : 1 : 1 one finds that the above condition 
holds when a ^ 0.45, in other words, when chain gap is 
at least half the maximal gap. As discussed above, a 
reasonable estimate of the ratio a is close to one half, so 
it looks like the condition for existing of the 7r-shifts in 
grain-boundary junctions is barely satisfied. However, in 
this kind of experiments the effect must be very sensitive 
to the value of the chain gap. We will return to this issue 
later. 

Our final point concerns the experiments on twinned 
samples. Naively, one can assume that the OP in the 
chain bands have the same sign in all domains, and thus 
the tunneling currents from different domains cancel. It 
easy to see, however, that it is not the case. Supercon- 
ducting state in each domain is degenerate with respect 
to changing signs of the OP in all bands simultaneously. 
Thus the relative phase of the OP in neighboring domains 
will be set by the proximity effects. Arguments similar 
to those used above for the tunneling currents lead to the 
conclusion that OP in two adjacent domains with oppo- 
site orientations will have opposite signs of the OP in the 
same bands, thus maintaining the proper a/b asymmetry 
for the whole crystal. 

To summarize, we suggest that recent Josephson junc- 
tions experiments that discovered n phase shift between 
the tunneling currents in a and b directions in YBCO can 
be fully understood in terms of the pairing symme- 
try, when order parameters in bonding and antibonding 
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plane bands have opposite sign, provided that the chain 
band is properly taken into account. This model is able 
to explain non-zero tunneling current perpendicular to 
the planes, as well as independence of the experimental 
results on twinning. 

A "smoking gun" for this model would be an experi- 
ment on YBCO with the superconductivity in the chains 
intentionally destroyed by doping at the Cu sites (Fe, 
Ga), which should be compatible with the conventional 
s-pairing. A particularly interesting property of our 
"chain scenario" is that the experiments with the YBCO- 
conventional superconductor junctions should be 

substantially less sensitive to such doping than the exper- 
iments with the grain boundary junctions |^], since in the 
latter case the condition on the chain gap is much more 
severe. Interestingly, the only tunneling experiment on 
YBCO which indeed showed no 7r-shifts was that of d, 
which was using grain boundary junctions. This fact is 
a strong argument in favor of suggested model. 
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TABLE I. Partial contributions of the chain, plane- bond- 
ing and plane-antibonding bands to the density of states and 
plasma frequencies of YBa2Cu307 (from Ref. [M). 



bonding 

antibonding 

chain 



N(n) 
N 



23% 
55% 
22% 



60% 37% 
37% 15% 
3% 47% 



77% 
19% 

4% 



40% 
7% 
54% 



Fermi surface of YBa2Cu307, from Ref. 
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